LINEAR  FEATURE  DETECTION  IN  SAR  IMAGES 


Ph.  Courmontagne 


ISEM 


Laboratoire  L2MP  UMR  6137  CNRS 
Maison  des  Technologies 
Place  Georges  Pompidou 
83000  TOULON 
FRANCE 

philippe . counnontagne@isem.tvt . f r 


ABSTRACT 

The  problem  treated  in  this  paper  is  the  one  of  detection  and 
restoration  of  ship  wakes  in  Synthetic  Aperture  Radar  (SAR)  im¬ 
ages.  This  cannot  be  easily  made,  because  SAR  images  are  cor¬ 
rupted  by  a  granular  noise,  called  speckle  and  because  there  is 
no  information  about  the  direction  and  the  level  of  the  wake.  For 
these  reasons,  most  detection  algorithms  use  the  Radon  trans¬ 
form,  which  makes  square  a  straight  with  a  point  We  propose 
here  a  new  method,  based  on  the  manage  between  the  Radon 
transform  and  a  filtering  method  used  to  interpolate  the  image 
in  a  rotating  reference  system,  introduced  by  the  Radon  trans¬ 
form  theory.  This  filtering  technic  is  the  stochastic  matched  fil¬ 
tering  technic,  which  allows  to  maximize  the  signal  to  noise  ratio 
after  processing.  Experimental  results  on  SIR-C/X-SAR  images 
are  presented  and  compared  to  those  obtained  using  the  classical 
approach. 

1.  INTRODUCTION 

Recently,  a  great  deal  of  research  has  been  dedicated  to  ship 
wake  detection  in  Synthetic  Aperture  Radar  (SAR)  images.  In¬ 
deed,  it  is  well  known  that  SAR  images  are  able  to  show  ship 
wakes  as  lines  darker  or  sometimes  brighter  than  the  surround¬ 
ing  sea.  Most  of  the  detection  algorithms  use  the  Radon  transform 
[2,  7,  8].  Indeed,  when  an  image  contains  a  straight  line  or  a  seg¬ 
ment,  its  Radon  transform  exhibits  a  narrow  peak  if  the  line  is 
brighter  than  its  surroundings  and  a  trough  in  the  opposite  case. 
Thus,  the  problem  in  finding  lines  is  related  to  detect  these  peaks 
and  troughs  in  the  transform  domain.  Other  methods  use  the  detec¬ 
tion  on  both  ships  and  ship  wakes  [5].  Given  that  SAR  images  are 
affected  by  a  granular,  multiplicative  noise  (called  speckle),  most 
of  these  detection  algorithms  pre-filter  the  data  in  order  to  improve 
the  visibility  of  the  ship  wakes. 

We  know  that  the  application  of  the  Radon  transform  requires  the 
computation  of  interpolated  image  in  a  rotating  reference  system. 
In  this  paper,  we  propose  a  new  method  based  on  the  mariage  be¬ 
tween  the  Radon  transform  and  a  filtering  method.  We  use  this 
filtering  technic  to  compute  the  interpolations  of  the  SAR  image, 
in  order  to  estimate  properly  the  signal  of  interest  (the  ship  wake) 
in  the  rotating  reference  system.  This  processing  is  called  the 
stochastic  matched  filtering  technic  [1].  It  is  based  on  the  sig¬ 
nal  expansion  into  series  of  functions  with  uncorrelated  random 


variables  for  decomposition  coefficients.  This  corresponds  to  the 
Karhunen-Loeve  expansion  in  the  case  of  a  white  noise.  Because 
the  chosen  basis  functions  improve  the  signal  to  noise  ratio  after 
processing,  there  is  no  more  sinusoidal  curves  corresponding  to 
the  speckle  in  the  Radon  domain,  and  the  detection  of  the  peak  (or 
trough)  corresponding  to  the  ship  wake  is  improved. 

First  of  all,  we  recall  in  section  2  the  discrete  Radon  transform 
[4],  in  the  case  of  a  two-dimensional  signal.  Then,  we  present, 
in  section  3,  an  interpolation-filtering  method  for  noise  corrupted 
image  in  a  rotating  reference  system.  First,  we  recall  the  stochas¬ 
tic  matched  filtering  technic  and  then  we  describe  how  to  perform 
an  interpolation  based  on  this  method  and  using  the  discrete  co¬ 
sine  transform.  We  finish  this  section  with  the  explanation  of  the 
subimage  processing.  Next,  in  section  4,  we  propose  an  example 
of  application  of  our  processing  on  a  SIR-C/X-SAR  image,  which 
shows  a  moving  ship  and  its  dark  turbulent  wake.  We  finish  this 
article  with  a  comparison  of  our  results  with  those  obtained  by  the 
classical  approach,  which  uses  the  Radon  transform  based  on  the 
nearest  neighbor  interpolation. 

2.  THE  DISCRETE  RADON  TRANSFORM 

The  Radon  transform  on  Euclidean  space  was  first  established 
by  Johann  Radon  in  1917.  Nearly  half  a  century  after  Radon’s 
work,  the  Hough  transform  for  detecting  straight  lines  in  digital 
pictures  was  introduced.  But  this  transform  is  actually  a  special 
case  of  the  Radon  transform.  We  are  going  to  recall  in  this  section 
the  discrete  Radon  transform  of  a  two-dimensional  signal. 
Considering  an  image  /,  (M  +  1)  x  (M  -f  1)  pixels,  its  discrete 
Radon  transform,  1  is  expressed  by: 

Af/2 

I  {xe^B)  =  ^  I (xe  cos 9  —  ye  sin B, x$  sin B-\-ye  cos B) , 

where  xe  and  ye  are  integers  which  are  bounded  by  — ^  and 
B  corresponds  to  the  rotation  angle  and  takes  values  between  0  and 

TT. 

The  previous  equation  shows  that  the  computation  of  the  Radon 
transform  requires,  for  each  parameter  B,  the  calculation  of  the 
new  pixel  values  in  reference  system  Indeed,  the  new  coor¬ 
dinates  are  not  integer  values  and  so  are  not  corresponding  to  the 
native  mesh  of  the  image. 
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By  nature,  the  Radon  transform  accentuates  linear  features  in  an 
image  by  integrating  along  all  possible  lines.  The  result  is  that  an 
image  which  is  non  zero  in  a  single  point  (xo,  yo)  has  a  Radon 
transform  which  is  non  zero  along  a  sinusoidal  curve  of  equation 
xe  =  xo  cos  $  +  yo  sin  9.  The  phase  and  frequency  of  this  si¬ 
nusoidal  curve  depend  on  the  spatial  location  of  the  correspond¬ 
ing  point  in  the  original  image.  If  the  original  image  contains 
a  straight  line  or  a  segment,  its  Radon  transform  exhibits  a  nar¬ 
row  peak,  if  the  line  is  brighter  than  its  surrounding,  and  a  trough 
in  the  opposite  case.  The  coordinates  of  the  peak  (or  trough)  are 
(x^o ,  ^0 )  which  correspond  to  the  parameters  of  the  polar  equation 
xoq  =  X  cos  00 sin  Bo  of  the  straight  line.  Thus,  die  problem  of 
finding  lines  is  reduced  to  the  detection  of  peaks  and  troughs  in  the 
transform  domain.  The  Radon  transform  is  particularly  suited  for 
finding  lines  in  a  noise-corrupted  image,  because  the  integration 
process  tends  to  cancel  out  intensity  fluctuations  due  to  the  noise. 
For  this  reason,  we  can  find  in  the  literature  several  applications  in 
the  domain  of  wakes  detection  (see  [2, 7,  8]  for  example). 

3.  INIT:RP0LATI0N-FILTERING  OF  A  ROTATING 
IMAGE 

We  have  seen  that  the  Radon  transform  of  an  image  implies, 
for  each  value  of  the  9  parameter,  the  computation  of  an  interpo¬ 
lated  image.  We  want  to  use  the  Radon  transform  in  order  to  detect 
ship  wakes  in  SAR  images.  Given  that  SAR  images  are  corrupted 
by  a  granular  noise,  called  speckle,  it  is  of  great  interest  to  take 
into  account  this  noise  to  compute  the  interpolation  of  such  an  im¬ 
age,  in  order  to  give  a  good  estimation  of  the  signal  of  interest 
in  the  rotating  reference  system.  For  this  reason,  we  present  here 
an  interpolation  method,  based  on  the  stochastic  matched  filtering 
method,  which  pnnciple  is  to  expand  the  noise-corrupted  signal 
into  series  of  functions  with  uncorrelated  random  variables  for  de¬ 
composition  coefficients. 

3.1.  The  stochastic  matched  filtering  method 

Consider  a  two-dimensional  noise-corrupted  signal,  Z{x^y), 
defined  over  £)  =  |  -  T.  T]  x  |  -  T;  T].  This  one  corresponds  to  the 
superposition  of  a  signal  of  interest  S(x,  y)  with  a  noise  j5(x,  y): 

Z(x.  y)  =  S(  j.  y) -f  V(x,y)  G 

where  S(x.y)  and  B(x.y)  are  assumed  to  be  independent  and 
stationary. 

We  want  to  expand  simultaneously  the  signal  of  interest  and  the 
noise  into  senes  of  the  form; 


% 

5(x.y)=  Jim  Sn^n(x,y) 

n=  1 
.V 


B(x,y}=  hm  bn'^nix,y), 
.V-»oo  ' 
n=  1 


In  these  expressions,  ^n(x,y)  are  the  deterministic  linearly  in¬ 
dependent  basis  functions,  and  Sn  and  bn  represent  zero-mean, 
random  variables  expressed  by  the  following  relations: 


S{x,y)^n{x,y)dxdy 

B(x,y)$n(x,y)dxdy. 


(1) 


The  determination  of  these  random  variables  depends  on  the  choice 
of  the  set  of  deterministic  functions  {$n  (x,  y)}.  We  will  used  the 
set,  which  provides  the  uncorrelation  of  the  random  variables,  i.e.: 

r  {^nSm}  =  ■f' {Sn  } 

I  E{h^b,n}  =  E{bl)S„,„,. 

Now,  we  show  how  to  determine  these  functions  <^n(x,  y).  In  or¬ 
der  to  find  them,  let  us  consider  the  stochastic  matched  filtering 
technic  [1]. 

If  we  consider  a  deterministic,  stationary  two-dimensional  signal, 
called  S(x,  y),  which  is  defined  over  D,  cormpted  by  an  ergodic, 
stationary  noise  B(x,  y),  the  matched  filtering  technic  consists  in 
finding  a  function  $(x,y),  defined  over  D,  in  order  to  maximize 
the  signal  to  noise  ratio  K,  expressed  by  the  following  relation: 

^  {  !//d  r  } 

When  the  signal  is  not  deterministic,  but  a  random,  zero-mean, 
stationary,  two-dimensional  signal  5(x,y),  we  can  show  that  K 
can  be  explained  as  follows: 

{  i/Zd  r} 

{l/Zo  y)dxdyf'^ 

Given  that  this  signal  to  noise  ratio  can  be  rewritten  as  the  ratio 
of  two  quadratic  forms,  it  appears  to  be  a  Rayleigh  quotient,  so  it 
will  be  maximized  if  ^(x,  y)  is  the  two-dimensional  eigenfunc¬ 
tion  associated  to  the  maximal  eigenvalue  of  the  following  integral 
equation: 

//bTssCs  -  x',y  -  y')^n{x',y')dx'dy'  = 

(2) 

ffo  ^BB(.x-x',y-  y')^„{x' ,  y')dx'dy', 

for  all  (x,  y)  G  B  and  where  Tss  and  Tbb  represent  the  covari¬ 
ances  of  the  signal  and  of  the  noise,  respectively. 

Random  variables  Sn  and  bn  are  uncorrelated,  when  the  #n(x,y) 
functions  are  the  eigenfunctions  of  integral  equation  (2),  with  eigen¬ 
values  An  verifying: 

E{blY 


When  eigenfunctions  ^n(x,  y)  are  normalized  such  as  the  follow¬ 
ing  integral 


-  x',  y  -  y')$n(x,  y)^n(x  ,y')dxdxdydy 


takes  one  for  value,  we  can  show  that  functions  ^^n(x,  y)  are  ex¬ 
pressed  by: 


’®'n(x,y)  =  JJ  TBsix- x\y -y)^n(x\y)dx'dy. 


(3) 


In  these  conditions  and  considering  the  Zn  random  variables,  ob¬ 
tained  by  projecting  functions  $n(x,  y)  on  noise-corrupted  signal 
Z(x,  y),  we  can  show  that  the  use  of  the  following  expansion 


N 

Z(x,y)=  Um  Vz„4'„(x,j/) 

jV-^oo  ^ 


n=:l 
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corresponds  to  a  signal  to  noise  ratio  of  the  component  of 
Z(x,  y)  expressed  by: 


where  cr|  /cr^  is  the  signal  to  noise  ratio  before  processing. 

So  all  the  eigenfunctions  associated  to  eigenvalues  An 

greater  than  one  can  contribute  to  an  improvement  of  the  signal  to 
noise  ratio.  For  this  reason,  filtering  the  observed  signal  can  be 
made  by  keeping  all  the  components  with  a  signal  to  noise  ratio 
greater  to  a  certain  level,  anyhow  greater  than  1. 


3^.  Interpolation  using  the  stochastic  matched  filtering  method 

To  compute  an  interpolation  based  on  the  stochastic  matched 
filtering  method,  the  basic  idea  is  to  expand  the  observed  signal 
and  then  to  restore  the  signal  of  interest  using  the  y)  func¬ 
tions  previously  interpolated.  But  this  reasoning  presents  some 
defaults.  Indeed,  it  implies  an  heavy  CPU  budget  and  some  mem¬ 
ory  problems  may  appear.  For  these  reasons,  we  are  going  in  this 
section  to  propose  a  new  formulation  for  the  stochastic  matched 
filtering  method  by  using  the  discrete  cosine  transform. 

3.2.  L  Analytical  approximation  for  the  solutions  of  the  inte¬ 
gral  equation 

We  can  find  in  the  literature  several  works  based  on  the  stochas¬ 
tic  matched  filtering  technic  in  its  discrete  form  (see  [6]  for  exam¬ 
ple),  Unfortunately,  when  we  consider  the  discrete  form  of  in¬ 
tegral  equation  (2),  the  eigenvectors  solution  of  this  generalized 
eigenvalue  problem  are  linked  to  the  native  increment  of  the  im¬ 
age  and  could  cause  problems  for  image  interpolation.  So,  we  do 
not  consider  the  discrete  relation  but  the  continuous  one  to  find  the 
$n(x,  y)  functions. 

Considering  the  DCT’  coefficients 

functions  ^n{x,y)Sss{x-x\y-y')^nATBB{x-x%y-y), 
we  can  show  that  solving  integral  equation  (2)  becomes  equivalent 
to  solving  the  following  linear  system: 


3.2.2.  Interpolation-filtering  method 

We  propose  in  this  subsection  a  new  formulation  of  the  stochas¬ 
tic  matched  filtering  method  by  using  the  discrete  cosine  trans¬ 
form.  In  these  conditions,  we  are  looking  for  the  expression  of  the 
coefficients  of  the  filtered  signal  expanded  into  cosine  series;  we 
shall  reconstruct  the  approximation  of  the  restored  signal  in  the  fi¬ 
nal  phase  of  the  processing. 

Let  Z(x,  2/)  be  the  observed  signal  to  be  expanded  and  let  Z(x,  y) 
be  the  reconstructed  filtered  signal.  We  have: 

Q 

z(x,  y)=^y2  y)  y)  ^  r>,  (S) 

n=l- 

where  Q  is  chosen  such  as  Xq  is  greater  to  a  certain  threshold, 
anyhow  greater  than  1. 

In  the  last  relation,  Zn  are  the  random  variables  to  be  determined 
from  the  input  data: 

~  if  (6) 

We  have  seen  that  these  random  variables  are  uncorrelated  when 
functions  ^n(z:,  y)  are  the  eigenfunctions  of  integral  equation  (2). 
The  ^n{x,  y)  basis  functions  are  obtained  by  projecting  functions 
^n{x,  y)  on  the  noise  covariance  as  described  in  relation  (3). 

First,  we  modify  this  relation,  in  order  to  express  the  coeffi¬ 
cients  of  the  ^n(rr,  2/)  cosine  series.  It  comes: 

Nf  Nf 

(7) 

it=0  1=0 

foT(p,q)  =0,  1,  Nf. 

In  like  manner,  from  expansion  (5),  we  obtain  for  expression  of 
the  DCT  coefficients  of  restored  signal  Z(x,  y): 

Q 

Sk,i=y2zr,^i,.  (8) 

n.=:l 


Nf  Nf  Nf  Nf 

E  E  = An  E  E 


(4) 


To  end,  we  have  to  explain  the  Zn  coefficients  in  terms  of  the  co¬ 
efficients  of  the  observed  signal  and  of  the  eigenfunctions.  We  can 
show  that  relation  (6)  is  equal  to: 


where  {Nf  -f  1)  corresponds  to  the  number  of  DCT  coefficients 
taken  account  and  is  high  enough  to  ensure  the  uniform  conver¬ 
gence  of  the  series  to  their  respective  functions. 

Finally,  the  analytical  approximation  $n(:r,2/)  of  the  ^n{x,y) 
functions  solution  of  the  integral  equation,  is  obtained  with  the 
following  relation,  for  (x^y)  €  D: 


This  new  method  for  finding  an  analytical  approximation  for  the 
solutions  of  the  integral  equation  has  been  quantified  and  com¬ 
pared  to  the  classical  approach,  in  the  case  of  the  Fredholm  inte¬ 
gral  equation,  that  is  when  the  noise  covariance  describes  a  white 
noise  [3]. 

^  Discrete  Cosine  Transform 


Nf  Nf 

<9) 

p=0  q—O 

where  represent  the  Z(x^  y)  DCT  coefficients. 

Considering  the  0  rotation  angle,  it  is  possible  to  compute  the 
interpolated-filtered  signal  using  the  following  relation: 

rr  /  \  V^V^'q  fTrk{xe-^T)\  f7cl{yd—T) 

Zeixe,  ye)  =  2^2^  (  oT - -  )  V  "2? — ~ 

where  {xe,ye)  represents  the  new  coordinates  of  the  pixels  in  ro¬ 
tating  reference  system 

With  such  a  formulation,  it  is  possible  to  use  the  stochastic  matched 
filtering  method  for  image  interpolation  in  a  very  short  time,  be¬ 
cause  all  the  computations  can  be  made  by  using  only  the  algo¬ 
rithm  of  fast  discrete  cosine  transform. 
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33.  Subimages  processing 


interest  5(x,  y).  We  can  show: 


To  apply  the  stochastic  matched  filtering  method,  it  is  neces¬ 
sary  to  respect  the  stationary  condition  for  the  signal  and  the  noise. 
We  know  that  we  cannot  consider  an  image  as  a  realization  of  a 
stationary  process.  But  after  segmentation,  we  can  define  several 
areas  representative  of  a  texture.  So  a  particular  area  of  the  im¬ 
age  (a  subimage)  can  be  considered  as  a  stationary  process.  For 
this  reason,  we  are  going  to  cover  the  noise-corrupted  image  with 
subimages.  The  subimage  size  is  chosen  such  as  the  subimages  are 
assumed  to  be  a  texture.  For  each  angle  we  apply  the  proposed 
processing  on  subimages,  with  M  x  M  pixels.  Scanning  ^1  the 
image  allows  a  complete  processing. 


Figure  1 :  Subimages  processing  showing  the  overlapping  between 
adjacent  subimages 


When  we  compute  the  interpolation  of  a  subimage,  we  find  an  edge 
effect  Indeed,  in  rotating  reference  system  the  coordinates  of 
the  pixels  localized  on  the  edge  of  the  subimage  do  not  depend  on 
the  native  subimage  in  reference  system  This  edge  effect  being 
maximal  for  angle  $  equal  to  the  subimage  size  after  process¬ 
ing  is  X  iV,  with  N  equal  to  To  limit  these  edge  effects, 
it  is  necessary  to  segment  the  native  image  to  obtain  subimages 
which  overlap,  as  shown  in  figure  1.  The  gray  areas  correspond  to 
the  superposition  of  subimages,  d  and  h  coefficients  represent  the 
distance  between  the  center  of  two  adjacent  subimages  (in  line  or 
row)  and  the  width  of  the  overtyping  area  respectively.  We  have: 

iVsin^ 
sin  0  -H  cos  ^ 

We  now  apply  the  interpolation-filtering  method,  proposed  in  the 
previous  section,  to  each  zero-mean  subimage.  Assuming  that  the 
noise  is  high-frequency  compared  to  the  signal,  when  we  apply 
this  processing  to  the  whole  noise-corrupted  image  with  the  same 
number  Q  of  basis  functions  for  each  subimage,  the  resulting  im¬ 
age  may  be  smoothed  or  still  noise-corrupted.  Indeed,  the  signal  to 
noise  ratio  is  not  the  same  for  each  area  of  the  image.  For  this  rea¬ 
son,  we  are  going  to  process  each  subimage,  with  different  number 
Q  of  basis  functions.  To  find  this  number,  let  us  consider  mean 
square  error  e  between  reconstructed  signal  Z(x,  y)  and  signal  of 


d  = 


N 


cos  0  +  sin  0 


and  h  ’ 


Q  Nf  Nf 

€  =  <Ts  +  j  ^  (<Tb  -  Ancrl)  ^  ^  • 

71=1  fc=o  i=o 

So,  for  each  subimage,  we  compute  e  for  different  values  of  param¬ 
eter  Q  (Q  being  in  the  interval  [1;  Qmax]>  such  as  greater 

than  1).  We  only  keep  the  Q  walues,  which  minimize  the  mean 
square  error.  The  noise  power,  0% ,  will  be  computed  in  a  homo¬ 
geneous  area  of  the  whole  noise-comipted  image  and  the  signal 
power,  0-5,  will  be  estimated  in  the  subimage  to  be  processed. 
Considering  now  the  problem  posed  by  the  overlapping  areas,  the 
corresponding  processed  pixels  will  be  computed  by  averaging 
each  pixel  having  the  same  position  in  the  overlapping  subimage. 

4.  SHIP  WAKES  DETECTION 

To  illustrate  our  processing,  we  have  chosen  to  apply  it  to 
an  image  acquired  by  the  Spacebome  Imaging  Radar-C/X-band 
Synthetic  Aperture  Radar  (SIR-C/X-SAR),  which  shows  a  moving 
ship  and  its  dark  turbulent  wake.  This  image  is  presented  figure  2. 


Figure  2:  SIR-C/X-SAR  image  (698  x  698  pixels) 


The  image  size  is  698  x  698  pixels.  The  number  of  gray  levels  is 
256  (0:  black,  255:  white).  The  dark  patches  in  the  upper  right 
of  this  image  correspond  to  smooth  areas  of  low  wind.  The  ship’s 
wake  is  about  28  kilometers  (17  miles)  long  in  this  image  and  in¬ 
vestigators  believe  that  may  reveal  that  the  ship  is  discharging  oil. 
Classically,  to  quantify  the  perturbation  level  of  a  SAR  image,  we 
determine  its  speckle  level.  This  one  is  obtained  by  computing 
the  variation  coefficient  (C  in  the  following),  obtained  on  several 
homogeneous  areas  of  the  image.  Let  W  be  the  number  of  homo¬ 
geneous  areas  /n,  we  have: 


1 

_ 
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For  the  studied  image,  the  variation  coefficient  is  equal  to  0.277. 
We  have  now  enough  information  about  the  studied  image  to  pro¬ 
cess  it 


4.1.  Signal  and  noise  auto-correlation  functions 

We  have  seen  that  the  interpolation-filtering  method,  presented 
in  section  3,  requires  the  a  priori  knowledge  of  the  signal  and  the 
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noise  auto-correlation  functions.  In  order  not  to  favor  any  particu¬ 
lar  wake  orientation,  we  have  chosen  to  represent  the  signal  auto¬ 
correlation  function  by  an  isotropic  model.  This  one  results,  for 
different  B  values,  from  the  averaging  of  several  auto-correlation 
function  computations  of  a  two  colors  straight  line  placed  in  a  ro¬ 
tating  reference  system.  This  model  is  presented  figure  3.a. 


Figure  3:  Normalized  model  for  the  signal  and  noise  auto¬ 
correlation  functions  (T  =  1) 

The  noise  auto-correlation  function  is  presented  in  figure  3.b.  This 
model  has  been  obtained  by  averaging  several  realizations  of  noise 
auto-correlation  functions  computed  in  some  homogeneous  areas 
of  the  native  image. 

4.2.  Radon  domain  and  wakes  restoration 

After  zero-meaning  the  image  presented  figure  2,  we  have  pro¬ 
cessed  it  with  a  subimage  size  equal  to  17  x  17  pixels,  to  respect 
the  coherence  length  of  the  noise. 

We  present  in  figure  4  the  interpolated-filtered  image  for  angle  9 
equal  to  35® .  For  this  image,  number  Q  of  basis  functions  is  in¬ 
cluded  between  1  and  13  depending  on  the  native  signal  to  noise 
ratio  of  the  subimage  to  be  processed  (near  13  basis  functions  for 
the  restauration  of  the  wake  and  1  basis  function  for  the  rest  of  the 
image). 
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Figure  4:  SIR-C/X-SAR  image  in  reference  system  §^35,  (973  x 
973  pixels) 

Analyzing  this  figure,  we  see  that  the  proposed  processing  allows 


a  great  reduction  of  the  speckle  and  a  good  restitution  of  the  signal 
of  interest  (the  wake).  Indeed,  there  is  no  more  dark  patches  and 
the  variation  coefficient  is  now  equal  to  0.016,  so  there  is  an  im¬ 
provement  by  a  factor  of  18  of  the  speckle  level. 

We  present  in  figure  5  the  resulted  image  in  the  Radon  domain. 
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Figure  5:  Radon  domain  obtained  using  the  interpolation-filtering 
method 

In  the  transform  domain,  the  vertical  axis  represents  the  orienta¬ 
tion  of  each  integration  line,  while  the  horizontal  axis  represents 
the  distance  of  each  line  from  the  center  of  the  image.  The  trough 
corresponding  to  the  wake  is  clearly  evident  Its  vertical  position 
is  near  51®  and  corresponds  to  the  orientation  of  the  wake.  Fur¬ 
thermore,  several  sinusoidal  curves  of  poor  amplitude  regarding  to 
the  trough  amplitude  are  visible  in  the  Radon  domain  and  corre¬ 
spond  to  the  Radon  transform  of  the  residual  perturbations  after 
interpolation-filtering. 

From  this  transform  domain,  we  have  used  the  inverse  Radon  trans¬ 
form  to  find  the  location  of  the  wake  in  the  spatial  domain.  This 
inverse  transform  has  been  applied  to  the  image,  presented  in  fig¬ 
ure  5,  before  raised  to  the  power  of  three  in  order  to  improve  the 
amplitude  of  the  trough  in  regards  to  the  rest  of  the  image.  The 
different  interpolations  have  been  made  using  a  nearest  neighbor 
interpolation.  We  present  in  figure  6  the  resulted  image. 


Figure  6:  Restored  wake  from  figure  5  (696  x  696  pixels) 

The  resulted  image  shows  that  our  processing  allows  a  great  im- 
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provement  of  the  wake  readability.  All  the  disruptive  pixels  have 
disappeared.  The  variation  coefficient  for  this  image  is  equal  to 
0.002  compared  to  0.277  for  the  original  image. 

5.  COMPARISON  WITH  CLASSICAL  PROCESSING 

We  present  figure  7  the  transfonn  domain  of  the  SIR-C/X- 
SAR  image.  Each  image  in  the  rotating  reference  system  has  been 
computed  using  a  nearest  neighbor  interpolation. 


Figure  7:  Radon  domain  obtained  using  the  classical  approach 

In  the  transform  domain,  we  can  see  several  overlapping  sinusoidal 
curves.  They  are  due  to  the  presence  of  the  speckle  in  the  native 
image  and  they  limit  the  visibility  of  the  trough  corresponding  to 
the  wake.  Furthermore,  the  dark  patches  in  the  original  image  cre¬ 
ate  in  the  transform  domain  a  large  sinusoidal  curve  with  an  am¬ 
plitude  higher  than  the  one  of  the  trough.  So,  using  a  threshold  to 
extract  the  trough  will  be  misleading,  because  the  minima  of  the 
transform  domain  do  not  correspond  to  the  trough.  In  this  case,  the 
detection  cannot  be  released  with  this  processing. 


.4Kn  .too  e  3M 


Figure  8:  Restored  wake  from  figure  7  (698  x  698  pixels) 

To  illustrate  the  previous  remarks,  we  present  in  figure  8  the  result 
obtained  after  having  applied  the  inverse  Radon  transform,  with 
the  same  process  as  for  the  image  presented  in  figure  6.  It  confirms 
that  this  processing  is  not  efficient  for  such  an  image.  Indeed,  we 


remark  a  lot  of  perturbations,  mainly  due  to  the  presence  of  the 
dark  patches.  The  variation  coefficient  for  this  image  is  equal  to 
0.145  and  the  improvement  is  only  a  factor  of  2. 

6.  CONCLUSIONS 

We  have  presented  in  this  paper  a  new  processing  which  allows 
ship  wakes  detection  in  SAR  images.  This  processing  is  based  on 
the  computation  of  the  SAR  image  Radon  transform.  The  original 
contribution  of  this  work,  compared  to  the  classical  approaches 
in  this  domain  consists  in  taking  into  account  the  noise  for  the 
image  interpolation  in  the  rotating  reference  system.  This  allows 
the  perturbations  to  have  a  lower  impact  in  the  transform  domain, 
the  corresponding  sinusoidal  curves  having  an  amplitude  smaller 
than  the  p^  or  trough  characteristic  of  the  wake.  We  have  com¬ 
pared  the  transform  domain  and  the  restored  wake  obtained  by  our 
processing  on  SAR  images,  with  those  obtained  with  the  classical 
processing.  In  all  cases,  our  processing  presents  far  better  results. 
With  our  processing,  the  probability  of  false  alarm  or  no  detection 
is  lower  than  with  the  classical  approach,  because  only  the  signal 
of  interest  is  considered. 

An  important  drawback  of  the  Radon  transform  is  that  it  is  global 
by  nature,  so  this  transform  carmot  tell  the  difference  between  long 
and  short  straight  lines.  For  this  reason,  future  work  in  this  domain 
concerns  the  application  of  the  proposed  interpolation  method  for 
the  computation  of  the  localized  Radon  transform  [2],  which  al¬ 
lows  to  localize  the  beginning  of  the  wake. 
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